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ABSTRACT 


The reorder point and reorder quantity for a multi- 
echelon inventory system consisting of two levels were 
determined through the use of a mathematical model anda 
computer simulation was used to verify the results. 

The main echelon supported two lower echelon stock 


points and reordered using a continuous review inventory 


policy. The two lower echelon stock points operated under 


periodic review policies. 

The measure of effectiveness used was to minimize total 
system costs subject to a constraint on the maximum number 
of back orders per year. 

The results from the mathematical model were used as 
input values for the simulation and measures of effective- 
ness were compared. An alternative procedure was proposed 
and simulated; and the results of the three products were 


compared. 
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I. INTRODUCTION 


Despite the fact that all military supply systems and 
many large civilian corporations have multi-echelon and 
multi-location inventory systems, few multi-echelon inven- 
tory models are currently being used. The policies that 
have been used are the single echelon and single location 
inventory policies that attempt to minimize the total 
variable cost of a single location ignoring total system 
cost. The difficulties of optimizing a multi-echelon and 
multi-location inventory system are most likely due to the 
complexity of demand and the interdependency among the units 
in different echelons. 

A mathematical model for a multi-echelon and multi- 
location inventory system is developed in Chapter II. The 
objective of this model is to minimize the total annual 
variable costs subject to a maximum allowable number of 
back orders per year for the entire system. This model 
consists of one first echelon location operating under a 
continuous review policy and two second echelon systems 
which operate under periodic review policies. An example 
is presented to obtain numerical results. 

In Chapter III, a computer model is developed to simu- 


late the multi-echelon, multi-location system. This simu- 


lation model was used to check the results obtained from 


the mathematical model and to demonstrate the interaction 


between echelons. Graphical plots are generated which show 
the inventory position of each entity in the model. 

In Chapter IV, comparisons are made of the simulation 
results and the results obtained from the mathematical 
model. The comparisons suggest some modifications to the 
operating policy. 

The system was simulated again using the modified 
operating policy and substantial improvements in the 
effectiveness of the multi-echelon system were observed. 
These results are described in Chapter V. 


Chapter VI summarizes the results of the research and 


concludes with some suggestions for additional research. 


II. A PROBABILISTIC MATHEMATICAL MODEL FOR MULTI-ECHELON 

INVENTORY SYSTEM 
A. DESCRIPTION OF MODEL 

In this model there are three systems. One of these 
systems represents the highest echelon and is called the 
Main System. The other two, at separate locations, are 
the lower echelon and are called System One and System Two, 
respectively. Each system carries its own inventory, and 
receives random demands. System One and System Two are 
dependent upon the Main System but independent of each 
other. System One and System Two can only be resupplied 
by the Main System. In other words, they cannot order 
from any other suppliers. The Main System replenishes 


stocks by placing orders to external suppliers. 


B. OBJECTIVE OF MODEL 

The objective of this model is to minimize the total 
expected annual variable costs of three systems subject 
to a specified expected number of total back orders per 
year. In fact, this is the same as minimization of total 
yearly variable costs subject to a specified minimum level 
of customer satisfaction. 

Based on these objectives, decision variables for each 
system will be calculated *o achieve optimum levels for 


the entire multi-echelon supply system. 
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C. ASSUMPTIONS OF THE PROBABILISTIC MATHEMATICAL MODEL 
1. Main System 

It is assumed that the Main System under considera- 
tion consists of a single installation which utilizes trans- 
action reporting. This system is governed by a (Q, r) type 
policy with back orders. 

The assumptions in addition to the continuous review 
assumption are: 

a. The cost of operating the information processing 
system is independent of Q (reorder quantity) and r (reorder 
point). 

b. The unit cost C of the item is a constant in- 
dependent of Q. 

c. The back-order cost is constant (I[), per unit 
back ordered regardless of the length of time the back- 
orders exist. 

d. There is never more than a single order out- 
standing. This assumption implies that when the reorder 
point is reached, there are no orders outstanding; there- 
fore, the inventory position is equal to the net inventory. 
Thus, the reorder point will be the same regardless of 
whether it is based on the inventory position or net 
inventory. 

e. Procurement lead times are independent and 
identically distributed (i.i.d) random variables with a 


gamma distribution. 


ll 


ee Se: Pe oot il bases ee 


£. All variables are treated as continuous. 


mean number of demands per year a constant Ant 


g. The demands are Poisson distributed with the | 


h. The reorder point, r, which is based on the 


inventory position is positive. 


With a back orders constraint, it is infeasible to 


wait until back orders exist before placing an order. 


Because of this and assumption (d) there will be no back 


orders outstanding at the reorder point. As was discussed 


in assumption (d), at the reorder point, the inventory 


position is equal to the on-hand inventory. 


For this model, any one of the three inventory 


cere parent > 


levels on hand, net, or inventory position can be used to 


define the reorder point; and the reorder point has the 


same value for any one of them. It should also be noted 


that to use the on-hand level it must be assumed that 


after an order arrives, it is sufficient to fill all back 


orders and raise the on-hand inventory level above the 


reorder point. If this ever failed to happen, the reorder 


point would never be reached again and the system would 


j continue to accumulate back orders indefinitely. 


When the reorder point is thought of in terms of 1 


{ 
{ the inventory position of the system, then assumption (d) 
i 


guarantees that the on-hand inventory will always exceed 


the reorder point when an order arrives; otherwise, it 


would. not be possible to have only a single order outstanding. 
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2. System One 


This system, which is at the lower echelon, is 
governed by a periodic review policy. The operating 
doctrine is the most widely used type of periodic review 
which is the order up to R policy. All demands which 
occur when the system is out of stock are back ordered. 


For this periodic review system the time between 


reviews will be denoted by T, and at each review time a 
sufficient quantity is ordered to bring the inventory 
position of the system up to a level, R, regardless of the 
amount of on-hand inventory. This policy dictates that 

at review times even if the inventory position is R-l, 
only one item must be ordered to bring the inventory 
position up to R, ignoring the high cost of placing the 
order. Despite its appearance of being illogical, this 


assumption simplifies the formulation of the system and 


is very unlikely to occur on high demand items. When an 
item experiences zero demand in a review period, there is 
no need for order because the inventory position is 
already at R. 

The other assumptions are: 


a. The cost, J, of making a review is independent 


of the variables R and T. 
b. The unit cost, C, of the item is constant and 


independent of the quantity ordered. 
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c. Back orders are incurred only in very small 


quantities. This implies that when an order arrives, it 
is almost always sufficient to meet any gutstanaing back 
orders. 

dad. The backorder cost is constant, Il, per unit 
back ordered regardless of the length of time the back 
order exists. 

e. Procurement lead times are i.i.d random 
variables with a gamma distribution. 

f. Orders are received in the same sequence in 
which they were placed. It should be noted that for (Q, r) 
models, the two assumptions that orders were received in 
the sequence placed and that lead times are i.i.d random 
variables could not both hold rigorously, since there 
exists a positive probability that two successive orders 
could be separated by an arbitrarily short time interval. 
In this model, orders can never be more closely spaced 
than by an interval of length T. If T is large enough, 
it is possible, provided that there is a sufficiently 
small range of variation in the lead time, that both 
assumptions hold simulténeously. 

g- The demands are Poisson distributed with 


mean, A. 


h. All variables are treated as continuous. 
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3. System Two 


This system is identical to, but independent of, 


System One. 


D. ANALYTICAL SOLUTION TO THE MODEL 

As was mentioned previously, the objective of the model 
is to minimize total system costs subject to a constraint 
on the number of back orders per year. 

The determination of various annual costs of each 
system may be made independently and placed into a common 
cost formula. Instead of doing this independently, it is 
better to find a main system cost expression, System One 
cost expression and System Two cost expression, and then 
to add them to each other to determine the total cost 
formulation and then use this formulation as the objective 
function of total system. The systems are actually tied 
together through the constraint on the number of back 
orders. 

For the following notation, the subscript m will 
indicate that the variables being subscripted belong to 
the Main Sytem where 1 and 2 indicate System One and 
System Two, respectively. ! 

The costs that are of interest in this model for each 
system are the cost of placing an order, the cost of 


carrying inventory and back order costs. 
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1. Main System 


As was discussed previously, the Main System has 


a transactions reporting policy or (Q, r) model. 


In a continuous review system, a period is defined 
as the length of time between the receipt of two successive 
procurements. This time period is a random variable. 
Because the procurement lead times are random variables, 
the number of demands for a fixed time are random variables; 

j and the number of items demanded per demand are also random 
variables. The review period is also a random variable. 
In the following material, we itemize the costs. 


a. Procurement Cost 


: mn = Expected number of demands per year 
} 
i Q, = Order quantity 
rin 
on Average number of procurements per year 
m 
An = Procurement cost per cycle 
An 
a An = Procurement cost per year 
™ 


b. Inventory Carrying Costs 
Because of randomness, it is possible to accumu- 
late a large number of back orders at the end of a cycle. 
To prevent this from occurring, it is advisable to provide 


a safety stock to buffer the system from excessive numbers 
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of back orders. The safety stock is the expected amount 
of stock on hand when an order arrives. The actual amount 


of stock on hand when a shipment arrives is clearly random. 


S_ = Mean value of on-hand stock when an order 
arrives 
After an order arrives, the expected on-hand inventory in- 
creases to Q+S and is reduced to a value of S on the average 
just before the next order arrives. Therefore, the average 
on-hand inventory per cycle is: 


(Q,t8,) Sa Q 
Pom goo? © aa 2a 


To write Sin in terms of the reorder point, r., 
let us first assume that the lead time Th is fixed. Let 
& (xix) =r, 7 * be the net inventory 


at the time an order arrives, 
where x 


is the number of units demanded in lead time Tat Then, 


S_=E é {Net inventory] = Bam (Xi Em) £(xit.) dx 


where £(x;T)) = Density function of demand in time con 


Ss Beh (r-%) £(x;T) dx = In a f(xit) ax 


S. =r -p where a Expected lead time demand. 


Then = en =e FE 7 Um 


H 
i) 


Inventory carrying charge 
C = Cost of an item 


. _ = Qm 
Total holding cost/year = IC ea ee =) 


c. Stockout Costs 1 


0 if x-r<0 
Let us define n. (x;r)) = { 
m x-r Lf x-r 20 


where n. (x ;t)) is the number of back orders per cycle. 
m 


If An (Xm) = expected number of back orders per cycle, then 


Ma(Fm) 7 nN, (*;r,,) A(x) dx = ee (x-r) h(x) dx 


m 
ae xh (x) dx-r H(r,) where H(r)) = P[xX>r] 
m 
and h(x) = marginal distribution of leadtime demand. 


Therefore, the expected number of back orders/year = 


Xm 
5 | 


om xh(x) dx ~ xr H(r.)] 


co) 
x m 
m 


and the expected cost of back orders/year = 


1, 
mm 2 
—— lf xh (x) ax-r H(rJ)] 


om 


All the terms in the average annual variable 


cost Kn have now been found: 


r mA 
m Qn 


K_ = An + IC [>= + rma wal +—— [| a xh(x)dx-r_ H(r)] 


mo Qh m “m 


2. System One and System Two 


It was stated in the assumptions that System One and 
also System Two both follow a periodic review policy with an 
order up to R stockage policy. Since System One and System 
Two are identical to each other; the equations will be 
derived only for System One. 

For convenience, a period is assumed to be the time 
between the receipt of two successive orders rather than 
between the placement of two successive orders. Costs are 
described as follows: 

a. Ordering and Reviewing Costs 

J = Reviewing cost/cycle 


A = Ordering cost/cycle 
Ji +A 


T 


1 


Ordering and reviewing cost/year = 
1 


where Tj is the period length defined in units of years. 
b. Inventory Carrying Cost 
The expected net inventory just prior to the 


arrival of an order is Ry ns ae dy L’ where Mage mean demand 


during lead time. 
The mean rate of demand remains constant over 


time and the expected demand per period must be the expected 
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amount ordered, i.e., AT. If the expected net inventory 


immediately after the arrival of a procurement is Ry7 Hye 
it is therefore Ry - uy) > AT) just prior to the arrival 
of a procurement. 

The expected unit years of storage incurred 


per period is 


1 1 aay 
T, fg (Ry - wy) + F CR, - Wy - AQT))] = TLER] - by - —5>] 
BaP, 


and average inventory carrying cost/year = IC [R, - 


ec. Stockout Costs 

First we assume the case where the procurement 
lead time is constant <7. An order placed at time t will 
arrive in the system at time t + t, and the next procurement 
will arrive in the system at time t + Tt + Ti: After the 
order is placed at time t, the inventory position of the 
system is Ry: It is necessary to compute the expected 
number of back orders occurring between t + tT andt+%31+T 
A back order will occur in this period under assumption c if 
and only if the demand in the time period t + Ty exceeds R,- 
Assumption c also assures that after the arrival of the order 
placed at time t, there will be no remaining back orders, 


and therefore they must all occur between times t + 1 and 


t+t + Ti: Consequently the expected number of back orders 


incurred per period is 


zo (x-R, ) E(x; t + T,) ax where 


a 
Bsc ee T)) = Demand distribution during time Tt + Ty: 
When lead time is random’ with density g (4) 
with Tin, and Tas being lower and upper limits respectively 


and Tt, andT5, the lead times for the orders placed at times 
t and t + Ty respectively, the expected number of back 


orders incurred per period is: 


i f f (x-R, ) £ (x; T+T,)9(t4)g(t,) dx dt,dt, 


=f (X-R,) h(x; T,) ax 
1 


T 
max 
f 


where h(x; T)) = £(x; TotT, )g(t,) dt, which is the 


Tin 
demand distribution during time To + Ty when lead time is 
a random variable with density function g(t,). The average 


number of back orders incurred per year is; 


ss Ga oye n(x: 
Ej (R,,T) = t, Ri (Xx R,) h(X; T,) dx and the average 
back order cost per year equals toll, £, (R,T,) and 
n, (x) = Ty Ey (Ry To) is the expected number of back orders 


per period. 
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Finally, the annual variable cost of System One 


is: 
. AiT 


i 
Uy 2 


L 
1 
a Perth ] + 1, Ey (R,-T,) and 


likewise the annual variable cost of System Two is: 


L AaT. 


eS Shia ad 
Ko T + Ic [R, lg - Sy“) + I, EB, (R,,T,) 
where Li = J. +A.. 


3. Objective Function and Constraints 


The objective function of the model consists of the 


total annual variable costs of each system. Therefore; 


K = K, + kK, + K which is equal to the minimization 


1 2 
of; 
A Q Td L 
m ™ mm 1 
Kei A PEC: ee rd ie n (x) + = 
QO, Om 2 m m Oo, Um tT, 
A,T I L A5T 
2 ay Sal ee Som 2 ee | 
+ Ic (Ry Wy =] + Tt Ny(r) + Tt, + IC [Ro Wo" —5 
I 
2~ 
+ =~ ny (4) 
T, 2 
subject to im n(x) + a n, (cr) + = nj(r) <b 
J om ie soem le 
where b is the specified total maximum number of back orders af 


per year for the entire system. 

The problem at hand is to calculate the optimum 
values of Qua! Cn! Ri Tie Ro, To: Since the objective 
function and the constraint are non-linear functions of the 
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decision variables, we solve the problem using the Lagrange 
multiplier approach. 
4. Optimum Values of Operating Variables 
After including the Lagrange multiplier in the 


formulation, the new objective function becomes: 


Minimize 
. Q r 
om m mm - ye 
L==~A_ +Ic(-+= +r - ud + ce ie 
om m 2 m m Qn m Ty 
eT I L 44T 
Po - 11 es _2 4 _ 272) 
+ Ic(R, My 5 } + 7, ny (4) + T, + Ic{R, U5 5 
if r 
2 - m- 1 - 1 - 
+= ny (xr) - O((—— n(r) + S ny (r) + = Nn, (r)) - bd) 
TS 2 Qa m Ty 1 Ts Pa 


where 8 is the Lagrange multiplier. It is through this 
Lagrange multiplier that the three systems are linked together 
mathematically. 

The optimum values of the unknown variables can be 
found by taking the derivatives of the objective function 


with respect to Qn! Rae Ry, Tie Ro, Tos 8; equating them to 


zero; solving the equations simultaneously; and ensuring the 
Kuhn-Tucker conditions are satisfied. 


The derivatives were taken with respect to Ont Xm! 


Ry Ry 6. A different procedure was utilized to find T) 


and T5- 


The derivatives: 4 


A A 
GMs 21 MM Ic _ =m - a = 
™m m Q Q 
™m m 
L Ba eas Mie ge ceeds | 
oe 3. Qn m m 
m 
L 1 
L = —— = IC - = H,(R,,T,) (M1, - 8) = 0 
Ry Ry Ty hae eels L 
aL 1 
L =—— = IC - —H (R,, T.) (1, - 3) = 0 
Ry dR5 T, 2 2 2 2 
3L Ip i ge is 
Lo = se = (lm ni(r) + Sn (cr) + = n,(r)) - b] = 0 
8 38 Ci m T 1 T, 2 
The Kuhn-Tucker conditions: 
Lom 2 0 eg = 0 Qua 3 0 g <0 
Le. 2 0 in oe = 0 ee 0 §8.g = 0 
m m 
Lp 3 0 R,-Lp = 0 Ry 2 0 8 < QO 
1 1 
L 2 0 R,.L = 0 R, 3 0 
Ry 2 Ry 2 
x ‘ 
m= 1 - 1 -~ 
where g = =~ ni(r) + a n,(r) + = no(r) - b 
Qn m Ty ol T> 2 


Therefore the following Kuhn-Tucker conditions must 


be satisfied: 
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CER mea ss aN See ame ee ong 1 


~mm, Ic _ 
2 
Qn 2 Qn m 
ie 
TCs Bley) We 8) 80 
mi 
0 = Ae ays oh te ey Se 8 
se ea din” 29-2 
tc- + 8. (R., T.) (I,-8) 2 0 
eae aaa aaa 2 2 


m om m m= = _ 
oe: om + 2 om im (2) oes 8) 0 
Am m 
Ic r - H_ (r) (1. -8) = 0 
Qn 
Ry n 
Ic Ry a 7 Hy (Riy T,) (iy - 6) =0 
Ry ~ 
Qn 2 0 
rn > 0 
Ry > 0 | 
r] 
Ro > 0 
*m n (r) + ns n,(r) + ae (r) - b < 0 
Qn m Ty 1 Ts, 2 
‘nm = 1 ee 
8la- M(t) + He ny (¥) + a g(r) - b] = 0 
™ 1 
9 < 0 
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a idan acs 


ee ene nee ee 


From these equations and conditions the optimum 


values of the unknown variables can be simplified to the 
following form and the exact values can be found by 


solving these equations simultaneously: 


ma Ic Qu 
H (r) = 
m C1 =8) An 

Qn = 24m An + "mn Tn - 3 
Ic 

. io, 

Ba Sie Byte ea: 

> Ic Ts 


The procedure to find T) and T5 and subsequently Ry 
and Ry is through iteration and trial and error. This is 
not unreasonable. In realistic cases other considerations 
not modelled here usually dictate the length of a period. 

In most cases the systems control a large number of different 
items and the same period length is used for each item. 

Thus, the period length is usually some convenient calendar 
or financial period, such as a month or a quarter. The 


procedure we use is shown in the following example. 


E. EXAMPLE PROBLEM 


l. Additional Assumptions 
For this example it is assumed that each system has 


Poisson arrivals independent of each other. The customers 
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demand only one item at a time. System One and System Two 


peer 8s 


are resupplied only from the main system and the main system 
is resupplied from outside suppliers. To be more realistic 
it is also assumed that the procurement lead times have a 


gamma distribution with different mean and variances. 


copes enero RD 


The A's will be assumed daily demand or arrival 

rate per day rather than yearly values. 
2. Input Values 

Before providing the input values, some clarifications 
must be provided regarding the demand at the main system. 
When it is stated that 1, = 4/day, Ay = 3/day and Ay = l/day. 
These represent the mean number of demands which arrive each 
day directly at the respective systems. However, since all 
demands at Systems One and Two eventually must filter up to 
the main system, the cumulative demand at the main system 
has an expected value of Ay + ho + din = 8. The main system 
supports not only Systems One and Two, but also has its own 
customer demands. The model assumes that the demand at the 
main system is the superposition of the direct demands at 


the main system and SyStems One and Two. However, since the 


actual demands placed on the main system in the lower echelon 


{ 
{ 
| 
are batched, the demand variability is much greater than i 

would be expected by the superposition process. We will 

evaluate the seriousness of our assumption about the demand 

process at the main system with the simulation model 

described in the next chapter. 
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In this example the input values for each system 
are; 
a. Main System 
A. = 8/day (considering the individual customer 
arrivals and the remainder being Systems 
One and Two average demands) 


c = $50/item 


I = 0.23 
Tin = $5 
An $500 


b. System One 


hy = 3/day 


c = $50/item 


I = 0.23 
y = $5 

Ay = $100 
Jy = $100 


c. System Two 


Ay = 1/day 


Cc = $50/item 
I = 0.23 
I, = $5 


Ay = $100 
J, = $100 
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The total system is expected not to exceed 
b = 175 total number of back orders at the end of a year. 
3. Solution Procedure and Results 
It was assumed that for each system, procurement 
lead times are gamma distributed random variables. 


The gamma distribution is defined as: 


t 
g(t) = on gent e 8 with 
T(a) B 


yw = a-B and 
Var = 1. B2 


We need the lead-time demand distribution for the 
case in which lead times are gamma distributed and the 
process generating demands is Poisson. This is derived 
below. 


The Poisson distribution is 


(xt) * et ; : ; 
£(x) = ane eS and the demand distribution 


during lead time is: 


£ (x;T) =f £(x) g(t) at 


_/ (at) * ett pond e t/8 se 
“9 x Tl{a) po 
x n -(F + At 
. tH Bax! ie ae "*e ‘ as 
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ee nN 
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The parameters « and 8 for each system are: 


u" 
N 
e 
aw 


ee 12.8 oy = 5.43 a 


a 3.125 By 


m 
il 


The main system has a Negative Binomial distributed 


lead time demand with parameters: 


8- (12.8) (3.125) = 320 


we 


Var = 8° (12.8) (3.125) [1+8-(3.125)] = 8320 


Cc 
i] 


System One has the parameters: 


= 3°(5.43)°(4.6) = 74.934 


var, = 3° (5.43) (4.6) [1+3(4.6)] = 1109.0232 


System Two has the parameters: 


= 1°(2.5)-6 = 15 


var, = 1° (2.5)-6 [1+1(6)] = 105 


Because the negative binomial is computationally 


intractable, we use the normal approximation for the calcu- 
lation of lead time demand probabilities. The normal a 
distributions are assumed to have the same mean and variance 
as the negative binomial distributions they replace. 

To solve this example problem, it is necessary to 


make an initial estimate for the Lagrange multiplier 4 
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which must be less than or equal to zero. After the initial 
estimate if the total number of back orders per year exceeds 
b, then the absolute value of 6 should be increased gradually 
until the number of back orders converges to b. 

For this problem 6= -1.5 works very well. For the 


main system the formulas are: 


Qn = mem * Myth) Chm ic 
| 
Ic 


O. is calculated to be 550 and 


= —_ 550 (0.23)°50 _ 
Hy f®) = Pp [XzrJ] = “TESL. 8) 3920 = 0.3332455216 


From the inverse standard normal distribution TI-59 
calculator program (Appendix G), the reorder level at the 


main system is found to be: 


r= dt (0- 4305333395) + = 359.2706827 


3 - 
f(r) = une) He) + o ee) 


where §@ is the functional value of standard normal distribu- 


tion at 0.4305333395 
Am *) = 20.08139334. Using this value in the G-. 


formula yields 
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Q. = 72.2920 {500 + 20.08139334 (5 + 1.5) 


™ = 568.84 


0.23). 
which is not equal to the first estimated value. If Qn is 


taken 569 then: 


369 (0-23) -50 _ 09, 3447576396 


H,(%) * TSF Ls) 3920 = 


ae (0.399071201) + He 356.4008941 


he (2) = 21.05458646 


Using this value in the Qn formula yields 


Qn = ¥2.2920 [500 + 21.05438646 (5 + 1.5) = 568.69 
(0.23) +50 


which is very close to the initial Qn value. A summary of 


the results for the main system follows: 


ot = 569 


S 36.400891 


m 


Am (x) = 21.05438646/period 
Tm (rt) = 108.047115/year 


Kn = $6796.50/year 


To solve for the optimum values for Systems One and 
Two a different procedure is followed. The total annual cost 
of Systems One and Two is a convex function of the period length 


T. For different values of T there are different total annual 
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cost values. The minimum of these values is the optimum 
total annual variable cost and the corresponding T and R 
values are the optimum operating values. A TI-59 program 
was written to perform the line search for the best value 
of T. The program is found in Appendix G. User information 
and the features of the program are also in the same appendix. 
This program evaluates the R value, safety stock, back 
orders per period, back-orders per year, annual reviewing 
and ordering cost, annual holding cost, annual back order 
cost, and finally, the total annual cost. 

Using this program the optimum values for Systems 


One and Two are: 


Ry = 306.809389 

T, = 2.39 months (72.6958 days) E 
S) = 13.78788901 | 
ny, (r) = 8.670289086/period 

Fy(r) = 43.53283223 

Ky = $2634.41/year 

Ro = 134.9807003 

T, = 4.07 months (123.7958333 days) sy 
S5 = 0 

7z(2) = 8.13480099/period 
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Ny (x) 23.98467122/year 


K $1377.56/year 


2 


The total back orders per year are: 


Nm (8) + ny (x) + No (x) = 108.047 + 43.532 + 23.984 


175.56 


which is the same as b = 175.56. 
The total annual variable cost is: 
K = Kn + Ky + K, = 6796.50 + 2634.41 + 1377.56 


K = $10808.47. 
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III. COMPUTER SIMULATION FOR MULTI-ECHELON MULTI-LOCATION 
SINGLE ITEM INVENTORY SYSTEM 


To check the analytical results, a computer simulation 
was written which uses the same operating assumptions and 
input parameters that were made for the analytical model. 
This model, however, simulates the real world more accu- 
rately since some of the simplifying assumptions required 
to obtain analytical results were not necessary in the 
simulation. Also, the demands placed at the main system 
from the lower echelon systems were batched as in the real 


world. 


A. DESCRIPTION OF THE SIMULATION MODEL 

As in the analytical case, there are three systems in 
the simulation model; the main system, System One and System 
Two. The flow charts of this program are in Appendix A. 

1. Main System 

This system uses a continuous review policy. When 

the stock level reaches the reorder point, it orders the 
amount Q. The decision variable for reordering is the in- 
ventory position. When an order is placed, the inventory 
position increases by an amount of Q. If the new inventory 


position is less than the reorder point, the system places 


an additional order for another amount of size Q. Then the 


inventory position increases with one more Q. This con- 
tinues until the inventory position exceeds r. The order 
policy is thus (nQ,r). N is the smallest integer that will 
make the inventory position higher than the reorder point. 

There are independent customer arrivals to the 
Main system and also group demands that are placed by 
System One and System Two when those systems replenished 
their stocks at the end of their periods. 

The number of units demanded per requisition is 
one. The number of units demanded for resupply to the lower 
echelon systems is random depending on the demands at the 
The 


lower echelon systems and the parameters R R 


1’ “2° 
program is capable of allowing geometrically distributed 
quantities demanded per requisition. Back orders at the main 
system are satisfied by filling the back orders to individual 
customers first followed by filling any back orders due to 
System One and System Two on a first-come, first-served 

basis. 

When a demand occurs from the lower echelon systems 
for which there is not sufficient on-hand inventory at the 
main system, the maximum amount is filled and the rest is 
put into the back order queue. As soon as a shipment 
arrives at the main system, the back orders are filled. 

The times between the customer arrivals are independent 
of each other and exponentially distributed. The lead 


times are also independent of each other and gamma distributed. 
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2. System One and System Two 
Both systems operate identically. The only differ- 


ence might be in the values of the system parameters. They 
have the same periodic review policy which is: order up to 
R at each review time. The program is also capable of making 
a decision at each review time regarding the placement of 
an order if the inventory position is less than or equal to 
a threshold value, is an (r,R) policy. If r is taken to be 
R-l, then the system orders up to R at each review time even 
if there is only one demand in a period. 

The times between customer arrivals are independent 
of each other and exponentially distributed. The number of 
units demanded per requisition is one. The lead times are 


independent of each other and gamma distributed. 


B. HOW TO USE THE PROGRAM 
Two Simulations are given in Appendix C and Appendix D. 
The programs differ primarily in the type of output that 
is generated. The program given in Appendix C produces a 
Versatec plot output showing the inventory position for 
the main system and the lower echelon systems. A simulated 
period of up to 4 or 5 years can be run with a Class K job. 
On the Versatec output for each system, there are two 
plots: one plot has a small triangle at each point repre- 
senting the inventory position. The second plot represents 


the net inventory. The first figure shows the main system 
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values, the second figure shows the System One values and 
the third figure shows System Two values all plotted versus 
time. 

The second program can be used to simulate operations 
over arbitrarily long periods of time. There are no dimen- 
sion restrictions. This program gives the net inventory 
at the end of each period for all the systems; total demand 
in the periods of both System One and System Two; demand 
during a lead time for the main system; average on-hand 
inventory; and average number of items short per day for 
each system. 

The input required by the programs is described by the 
variable definition list given in Appendix H. The programs 


use the IMSL subroutines for random number generation. 


Cc. RUNS 
1. Starting Conditions 

There were two different starting conditions entered 
for the main system. In the first case the inventory posi- 
tion and net inventory were set initially to equal the order 
quantity. In the second case, the inventory position and 
the net inventory were set equal to the reorder point plus 
the order quantity. For Systems One and Two the inventory 


positions and net inventories were set initially at R,- 
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2. Results of Runs 
In all runs the length of time was 10 years. The 
average yearly results were obtained by dividing the results 
for the 10 year period by 10. 
Using the policy parameters determined by the mathe- 
matical model and the first set of starting conditions the 


following simulated results were obtained. 


a. Main System 
The number of orders in 10 years = 5l. 
The number of back orders in 10 years = 3026. 
The average on-hand inventory over 10 years = 
313.026. 
The average safety stock when an order arrives = 
56.54 units. 


! The average annual costs are: 


Ordering costs: $2,550.00 
Holding costs: 3,599.80 
Back Order Costs: 1,513.00 


Total Cost $7,662.80 
b. System One 


Number .f back orders in 10 years = 587. 1 


The average on-hand inventory over 10 years 


114.482 units. 
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The average annual costs are: 


Ordering and Reviewing Costs: $1,004.18 


Holding Costs: 1,316.54 
Back Orders Costs: 293.50 
Total Costs $2,614.22 


c. System Two 
Number of back orders in 10 years = 287. 
The average on-hand inventory over 10 years = 
57.618 units. 
The average annual costs are: 


Ordering and Reviewing Costs: $ 589.68 


Holding Costs: 662.61 
Back Orders Cost: 143.50 
Total Costs $1,395.79 


(1) Total results 

The total number of back orders for the 
whole system per year is 302.6 + 58.7 + 28.7 = 390. 

The total annual variable cost for the 
whole system is $7,662.80 + $2,614.22 + $1,395.79 = $11,672.81. 

The model was run again with the second set 
of starting conditions and the same set of policy parameter 
values. The simulation results are summarized below: 

(a) Main System 

The number of orders in 10 years = 50. 


The number of back orders in 10 years = 


1665. 
41 


(b) 


(c) 


The average on-hand inventory over 


10 years = 346.937. 
The average safety stock when an order 
arrives = 77.46. 


The average annual costs are; 


Ordering Costs: $2,500.00 
Holding Costs: 3,989.78 
Back Order Costs: 832.50 

Total Costs: $7,322.28 


System One 

Number of back orders in 10 years = 582. 
The average on-hand inventory over 

10 years = 118.696 units. 

The average annual costs are; 


Ordering and Reviewing Costs: $1,004.18 


Holding Costs: 1,365.00 
Back Orders Costs: 291.00 
Total Costs: $2,660.18 


System Two 

Nunber of back orders in 10 years: 314. 

The average on-hand inventory over 10 years = 
57.401 units. 


The average annual costs are; 
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Ordering and Reviewing Costs: $ 589.68 


Holding Costs: 660.11 
Back Orders Costs: 157.00 
Total Costs: $1,406.79 


(2) Total Results 

The total number of back orders for whole system 
per year is 166.5 + 58.2 + 31.4 = 256. 

The total annual variable cost for whole system 
is $7,322.28 + $2,660.18 + $1,406.79 = $11,389.25. 

3. Comparison of Both Starting Conditions 
It is obvious that if both results are compared there is 

a considerable decrement on annual back orders but not a great 
difference in annual variable costs. 


390-256 _ | 
—<o0- 0.3436 


11672.81 ~ 11389.25 _ 
ema Co) 9 aaa ah 


43 


ne RO 


IN er 


Iv. COMPARISON OF ANALYTICAL AND COMPUTER SIMULATION RESULTS 


By comparing the analytical results with the simulation 
results, it is possible to evaluate the analytical model. 
If the measures of effectiveness predicted by the analytical 
model are reasonably close to those generated by the simula- 
tion, there is support for the analytical results. In the 
table below we compare the main measures of effectiveness: 


back orders and costs. 


Total Back Orders Total Variable 


Per Year Costs Per Year 
Analytical result 175.56 $10,808.17 
Simulation result 390.00 $11,672.81 or 
(For both starting 
conditions) 256.1 $11,389.25 


This table shows that the average yearly back orders are 
considerably higher than what is estimated by the mathematical 
solution. In fact, the total number of back orders generated 
in the simulation does not even meet the constraint. Since 
the solutions produced by the simulation model, especially 
with respect to the number of back orders do not agree with 
the mathematical model, this suggests that some of the 
assumptions made in deriving the mathematical results are not 
reasonable. The major differences between the simulation 


results and the analytical results are found in the measures 
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e 
of effectiveness at the main system. The results for the 


lower e®nelon track reasonably closely. 

The costs are also a little higher than what is 
- expected but much closer percentage wise. 

As an explanation for the large differences observed 
in the number of back orders at the main system with the ! 
simulation and analytical models, let us reexamine the 
assumptions we made that affect back orders. Back orders 
results from inadequate amounts of safety stock to protect 
the inventory system against excessively large numbers of 
demands in a lead time. The safety stock is manipulated 
by control of the reorder point r. After the reorder point 
is hit and an order placed, the system is totally at the 
mercy of the demands that occur during the lead time. If 
the variance of lead time demand is underestimated, then 
large stockouts will occur, even if the mean lead time 
demand is estimated accurately. The lead time demand is 
affected by two random quantities: (1) the demand distri- 
bution and (2) the lead time distribution. 

For purposes of making the mathematics tractable, we 
assumed in our analytical model that demands at the main 
system flowed in at a smooth continuous rate A which was 
taken to be the sum of the rates of demand incurred directly 


at the main system dan and those which occurred at the lower 


echelon systems Ay and Ag: The lead times were assumed to 


be normally distributed with mean equal to At and variance 

= An (1- 4g) where t is the mean value of the lead time. 
Because of the relatively high demand rates used in the 

example runs, the normal assumption should be justified by 

the Central Limit Theorem. However, as is illustrated by 

the Versatec plot output (See Appendix F) the demands at the 

Main system are far from smooth. What happens is the demands 

which arrive directly at the main system cause the inventory 

position to drop off smoothly. However, when the replenish- 

Ment orders from the lower echelon are received, large drops 

occur in the inventory position of the main system. Recall 

that the lower echelon systems order in batches once each 

period from the main system. If the main system simply tries 

to average out demands (as assumed by the analytical model), 

it will sometimes have very large amounts of excess stocks 

when shipments arrive and sometimes very large numbers of 

back orders. The high variance in lead time demand caused 

by the irregular demand actually seen at the main system 

causes the problem. Since all demands eventually flow 

through the main system and this is assumed by our analytical 

model, the problem cannot lie in the value used for the 

mean lead time demand. ' 


Let us explore further what happens when reorders are 


triggered at the main system. Demands directly at the main 


system eat away smoothly at the inventory position. Then a 
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very large quantity, say X, is demanded by one of the lower 


echelon systems. There is a very good chance that the 
large order placed by the lower echelon system will trigger 
a reorder by the main system. However, if the demand causes 
a large overshoot of the reorder point, the main system may 
have much less stock to live off of until the shipment 
arrives. For example, suppose the inventory position is 
IP = 347, the reorder point is r = 300, and System One places 
it's resupply order for 200 units. An order will be placed 
by the main system but instead of having 300 units to keep 
it going until the order is received, it will have only 
147 units. It is clear that if 300 is the amount of stock 
needed to provide reasonable protection against demands in 
a lead time, large numbers of back orders would be expected. 
Effectively, in the example above, the reorder level was not 
‘xr = 300 but r = 147 and the safety stock negative. The 
impact of this surge in demand caused by the batching of 
Cemands received from the lower echelon systems is to reduce 
the "effective" reorder point from the value r to a value 
r'<r. Because the actual demand distribution witnessed b, 
the main system is difficult to describe mathematically, 
the actual value of r' cannot be determined analytically. 
However, it is clearly less than r and may be much more so. 
In the next chapter we describe a policy modification 


for operation of the main system that was suggested by the 


observations above. The modification attempts to allow 


the main system to anticipate the surge of demands that 


will be received by the lower echelon systems. 
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V. AN ALTERNATIVE SOLUTION 


A. DESCRIPTION OF ALTERNATIVE PROCEDURE 

Due to the reasons mentioned in Chapter IV, the number 
of back orders found by simulation were much higher than 
the number of back orders predicted by the mathematical 
model. The question is how to run the multi-echelon system 
so that the large number of back orders seen earlier can 
be reduced. 

Obviously, what needs to be done is to reduce the impact 
of the very large demands that occur when the lower echelon 
systems place their resupply orders. With modern day 
communication and data systems, it would be feasible to 
allow the main system to "see" every demand that occurs 
anywhere in the system. If the main system is given the 
visibility, it will be able to take action to get the stock 
on its shelves in anticipation of the large demands from 
the lower echelon systems. The main system can do this if 
it uses as its reorder point the pseudo inventory position 
which is like the inventory position except that it decreases 
only when direct customer demands are encountered at any 


of the three systems. The pseudo inventory position is un- 


affected by the batch replenishment demands placed by the 


lower echelon systems. For example, if a customer requests 


a unit from System One, the pseudo inventory position at the 
main system decreases by one. However, if System One places 
an order for 200 units at the main system, the pseudo inven- 
tory position does not change. 

This new policy is referred to in this thesis as the 
"early warning policy." Clearly, the pseudo inventory posi- 
tion will always reach the reorder point before the inventory 
position. Therefore, orders will always be placed earlier 
and consequently, the number of back orders should decrease. 
The price paid will be in terms of extra holding costs. The 
results of a simulation using the early warning policy should 
be more nearly like those of the mathematical model since, 
in effect, the mathematical model makes the assumption of 
early warning. Implicitly, assumption of demand which is 
the superposition of the direct demands occurring at the 
Main system and Systems One and Two is equivalent to the 
"early warning assumptions." 

In the next section, results are given of a simulation 
of the multi-echelon system with early warning. The flow- 
chart of the simulation model is given in Appendix B. The 


actual FORTRAN computer program is given in Appendix E. 1 


B. COMPUTER SIMULATION RESULTS FOR EARLY WARNING POLICY 
The same starting conditions and input parameters used 
in the previous simulation run were utilized here. The 


results are summarized as follows: 
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l. Main System Results 


The number of orders in 10 years = 50. 
Total number of back orders in 10 years = 126. 
The average on-hand inventory = 488.093. 
Average safety stock = 210.96. 

The costs are: 


Ordering costs $2,500.00 
Holding costs 5,613.00 
Back order costs = 63.00 


Total costs $8,176.07 


2. System One Results 
Total number of back orders is 10 years = 665. 


it 


i The average on-hand inventory over 10 years = 117.846. 
j The costs are: 


Ordering and reviewing costs = $1,006.18 
Holding Costs = 1,355.23 
Backorder costs = 332.50 


Total costs $2,691.91 


3. System Two Results 
Total number of back orders in 1° years = 312. 
The average on-hand inventory over 10 years = 58.552. 
The costs are: 


Ope ee 


Ordering and reviewing costs = $ 589.68 
Holding costs = 673.35 
Back orders costs = 156.00 

Total costs $1,419.03 


C. COMPARISON OF EARLY WARNING POLICY SIMULATION RESULTS 
WITH ANALYTICAL AND FIRST SIMULATION RESULTS 


The comparison will be done with respect to total 
yearly back orders and total variable system costs (con- 


sidering all three system entities). 
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Entire System 


on a a 


Number of Back Total 
Orders Per Year Cost 
Analytical Result 175.56 $10,808.47 
First Simulation Results 256.1 11,389.25 
Early Warning Simulation 
Results 110.3 12,287.01 


Since the effects of implementing the early warning policy 
will be observed primarily at the main system, we also 


produce the results obtained for the main system individually. 


Main System 
Numbers of Back Total 


Orders Per Year cost 
Analytical Result 108.05 $6,796.50 
First Simulation Result 166.50 7,322.28 
Early Warning Simulation 
Results 12.6 8,176.07 


This table shows the differences better than the first one. 
The number of back orders decreases 88.3 percent, simultaneously 
as the costs increase about 10.44 percent. The reason for the 


higher cost is becausethe increment in safety stock increases 


the carrying costs. 
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VI. CONCLUSIONS 


The differences between the resultsof the mathematical 
and the simulation model can be explained largely as a result 
of the assumptions made about the demand process. In the 
main system it was assumed that the demand was smooth (the 
superposition of three Poisson processes), but in the simu- 
lation model the actual demand at the main system was as it 
would be in actual practice. There were batches of demand 
placed by System One and System Two in addition to the indi- 
vidual customer demands directly at the main system. These 
demand batches in fact increased the variance in lead time 
demand beyond that modelled. This explains why many more 
back orders were generated in the simulation model than 
what was predicted by the mathematical model. 

The simulation model developed in this thesis is useful 
for making comparisons and examining the effects of policy 
changes or parameter changes. Moreover, it is one of the 
best ways to check the reasonableness of all of the simpli- 
fying assumptions made in order to obtain analytical solutions. 

The mathematical results described in this thesis do 


not adequately determine the reorder point. The predicted 


number of stockouts is much less than the simulated numbers. 


As explained earlier, this is probably due to the assumptions 
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made in the model about the variance of lead time demand. 


Further study is needed to determine what could be done 
in the analytical model to better approximate what happens 
in actual practice. 

The early warning policy discussed in this thesis did 
provide for great reductions in the number of stockouts 


system wide. Since stockouts are probably the most important 


consideration in military supply systems, the early warning 1 


policy is recommended, even though the holding costs are 
larger. Additional study is required to see if the reorder 
level in the early warning policy can be reduced substantially 


from the value determined by the mathematical model. Pre- 


liminary evidence is that the reorder level can be reduced 


significantly (25 percent or so) without generating excessively 
many back orders if the early warning policy is used. 

Our simulation models allow us to view the effect of 
changes but they cannot be used to optimize the values of 
the policy parameters. For that objective, additional work 
in the mathematical modelling area is required. 

In this thesis we have tried to model a multi-echelon 
inventory system analytically, by linking together the 
individual echelons and locations through a single objective 
function and a constraint on backorders system wide. We, 
knowingly, were making various simplifying assumptions to 


facilitate the derivation of solutions. As reported above, 
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the resulting solution for the reorder level at the main 
system led to many more back orders than predicted. The 
other solutions; Ry: Ro, Q, The and T, appear to be satis- 
factory. 

In order to model adequately the multi-echelon system, 
it will probably be necessary to build into the determina~ | 
tion of the reorder lead at the main system the values of 
the parameters R and R, at the lower echelon systems. 

This will be the only way to accurately describe the actual 


demand process that is observed at the main system. We 


recommend future work in this area. 


APPENDIX A 


FLOWCHARTS OF FIRST SIMULATION MODEL 


Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 
Subroutine 


Subroutine 


Ship arrivals to System One 


One 
Two = Ship arrivals to System Two 
Three = Ship arrivals to Main System 


Four = Periodic review of System One 


Five = Periodic review of System Two 
Six = Shipment arrival to System One 
Seven = Shipment arrival to System Two 


Eight = Shipment arrival to Main System 
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Flowcharts of Early Warning Simulation Model 


Subroutine One Ship arrivals to System One 
Subroutine Two = Ship arrivals to System Two 


Subroutine Three = Ship arrivals to the Main System 


Subroutine Four = Periodic review of System One 
Subroutir.e Five = Periodic review of System Two 
Subroutine Six = Shipment arrival to System One 
Subroutine Seven = Shipment arrival to System Two 


Subroutine Eight = Shipment arrival to the Main System 
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APPENDIX C 


Simulation Program; Versatec Plotter 
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APPENDIX F 
Versatec Output of the Program in Appendix C 
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APPENDIX G i: 


TI-59 Calculator Program for Analytical Solutions for i 
Periodic Review Systems 


This program was written for probabilistic periodic 


review inventory models having gamma distributed lead times 
and Poisson arrivals and uses. Normal distribution for 
demand during lead time instead of Negative Binomial. 
Input Requirements 
The following variables should be stored in the regis- 
ters shown before the variables. 
STO 01 = A (arrival rate per day) 
STO 02 =C 
STO 03 = I 
STO 04 = I 
STO 05 = J 
STO 06 =A 
STO 07 = a-l 
STO 08 = 8 
This program also requires the TI-59 applied statistics 


module. i 
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Variable Definitionsfor Simulation Programs 


APPENDIX H 


Net inventory 

Inventory position 

Second inventory position for early warning system 
Order quantity 

Reorder level 

Number of items demanded by systems 

Number of items demanded per demand by ships 

Ship queue 


Lower echelon's group demand queue in the main 
system 


Amount of demand for each demand waiting in the 
queue 


Amount of shipment arrived 

Index. If it is equal to l, that means that the 
demand waiting in the main system queue to be 
filled belongs to System One. 

Total waiting time 

Total average on-hand inventory 

Length of a period 

This indicates the subroutines 

This indicates time of subroutines scheduled 
Net inventory variable for Versatec plotter 
Inventory position variable for Versatec plotter 


Time for Versatec plotter 
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Work space for geometric random variable 


Work space for geometric random variable 
Work space for geometric random variable 
Exponential random number 

Geometric random number 

Increment 


Indicates the number of ships waiting in the ship 
queue at System One 


Indicates the number of ships waiting in the ship 
queue at System Two 


Indicates the number of ships waiting in the ship 
queue at the main system 


Indicates the number of demand batches waiting 
in the group demand queue at the main system 


Time indicator 


Counter for ship arrivals to indicate the number 
of items demanded per period 


Multiplier for the number of batches of demand to 
be ordered from outside supplier 


Indicates the last change on net inventory of the 
main system in order to get the number of items 
demanded in a lead time 

Time of last change 

Scale parameter for a gamma distribution 

Shape parameter for a gamma distribution 


Ship arrival rate per day 


Probability of success for geometric distribution 
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